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Abstract 

In the said treatise Archimedes determines the equiUbrium positions of a floating 
paraboloid segment, but only in the case when the basis of the segment is either 
completely outside of the fluid or completely submerged. Here we give a mathe- 
matical model for the remaining case, i.e., two simple conditions which describe 
the equilibria in closed form. We provide tools for finding all equilibria in a reliable 
way and for the classification of these equilibria. This paper can be considered as 
a continuation of Rorres's article [8]. 

Introduction 

Archimedes' treatise "On Floating Bodies" (its customary Greek title is Uepl ^oxovfieuou, 
which literally means "about hovering things", see [T]) has been highly esteemed by 
mathematicians over centuries. Book 2 of this treatise can be considered as a sort of 
crown of Archimedes' work. In this book he applies a number of his principal results 
about volumes and centers of gravity to a problem which is extremely difficult to handle 
under Greek premises: namely, the determination of the possible equilibrium positions 
of a floating paraboloid segment (for details see Section 1 below). 

It seems, however, that even Archimedes was not in a position to treat this problem in 
full generality, since he restricts himself to the cases when the basis circle of the segment 
either lies outside of the fluid (i.e., the fluid touches this circle in at most one point) or 
is completely submerged. 

The case not considered by Archimedes occurs when the basis circle is partially sub- 
merged and partially not. It is, indeed, of a different nature than the "archimedean" 
case. Whereas Archimedes gives ruler and compass constructions for the "tilt angle" 
of the segment, results of this kind cannot be expected in the "non-archimedean" case. 
However, it is possible to establish a mathematical model for this case, based on two 
equations E = (the equilibrium condition) and F = (the floating condition). Unlike 
the corresponding equations in Archimedes' situation (see Section 1), E and F are no 
more purely algebraic expressions but also involve the function arctan(a;). Nevertheless, 
these expressions are rather simple if established with care (see Theorem [1]), and they do 
not involve "page-long monstrosities" of which Rorres [8] is warning. 

The said paper of Rorres contains a graphic completion of the "equilibrium sur- 
face" (based on numerical integration) together with interesting observations of physical 
phenomena and inspiring examples (one of which we repeat here in a treatment dif- 
ferent of his, see Example 1). However, [8J does not contain any closed formulas that 
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would describe the non-archimedean case (which is what we expect from a mathematical 
model). From this point of view j^8j appears only as a first step towards a completion of 
Archimedes' treatise. We hope that the present paper forms a second step. Such a step 
also requires simple tools by which one can decide whether an equilibrium position is 
stable or unstable. In this connection our above equations E = and F = are, again, 
quite helpful, since they imply that the Hesse matrix of the potential function looks fairly 
simple for equilibria (see Section 5). 

Maybe the most interesting question in the non-archimedean case concerns the number 
of possible equilibria for a paraboloid segment of a given shape and a given (relative) 
density. We have no mathematically rigorous answer to this question (which would 
constitute a third step). But we solve a related problem, namely, we determine the 
number of possible equilibria if the shape of the segment and the size of the submerged 
part of the basis are given, see Theorem O Combined with other devices like Proposition 
|2l this theorem allows finding "all" possible equilibria for a segment of given shape and 
density — not in the strict sense of the word but in a convincing manner, as we think. 

The standard English translation of Archimedes' treatise seems to be that of Heath 
[2] from 1897, which is based on Heiberg's first edition of the Greek text. The German 
version [3] (with helpful notes) is a translation of Heiberg's second edition [Ij, which 
considers the important Constantinople palimpsest discovered in 1899. We highly rec- 
ommend the monograph [3] about Archimedes and his works. For a survey of the treatise 
on floating bodies the reader may also consult [5]. A number of problems of floating ho- 
mogenous bodies are studied in |1] and f6j, works which also provide a wider theoretical 
background than we use here. For additional important references see [H]. 

1. The archimedean case 

Throughout this paper we denote subsets of in an abbreviated way; for instance 

{x < ay + b, y > c, z < dx^} 

stands for 

{{x, y, z) eM.'^ : x < ay + b, y > c, z < dx^}. 

In this section we give a modern paraphrase of Archimedes' results (see [1], Lib. II) and, 
thereby, introduce some basic notations. 

It suffices to consider the fixed paraboloid {z = x^ + y^} that arises from the parabola 
{y = 0, z = x"^} in the xz-plane by rotation around the 2;-axis. Our paraboloid segment 
V is defined by 

V = {x^ + y^ < z, z < a}, 

where a is the length of the axis {x = y = 0, < z < a} of the segment. 

Hence the parameter (in the usual sense) of the rotating parabola equals 1/2 and the 
geometric properties of V are completely determined by a. The surface of the fiuid is a 
plane 

£ = {z = bx + c} 

given by the parameters b, c. We assume b < throughout this paper. The segment V is 
said to be in right hand position ifV n{z > bx + c} lies outside the fluid (so in Diagram 



2 



1 this part of V is on the right hand side of £). Conversely, for a left hand position 
V P\{z <hx + c} must be outside. In order to exclude uninteresting cases one may also 
assume that both V ^{z > hx + c\ and V r\{z <hx + c} are non-empty. 





z 
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Diagram 1 



In the case considered by Archimedes the intersection of £ with the basis circle {x^ + 
< a, z = a} of the segment V consists of at most one point., which, in our context, 

must be the point (— \/a, 0, a). This situation is henceforth called the archimedean case. 

By our above "non-emptiness" assumption this case excludes a vertical plane £, so our 

definition of £ is sufficiently general. 

We treat the archimedean case for right hand positions first. To this end we note the 

volume V oiV and its center of gravity B, i.e., 

V = and B = (0, 0, 2a/3). (1) 

In order to define the axis of the paraboloid segment 

V' = Vn{z <bx + c} 

that forms the submerged part of V, we need the midpoint M = {b/2, 0, 6^/2 + c) of the 
line segment PQ, where {P, Q} is the intersection of £ with the parabola {y = 0, z = x"^}. 
A vertical line through M intersects this parabola in a point R (the vertex of V), and 
R = (6/2,0,6^/4). The axis of V is the line segment MR. One immediately finds that 
its length equals 

a' = by A + c. (2) 

Archimedes knew that V has the volume V = a'^7r/2 (see [T], Lib. II, Sect. iv). He also 
knew that its center of gravity B' = [x', 0, z') lies on MR in such a way that B'R has 
length 2a73 (ibid.. Sect, n); this gives x' = 6/2, z' = hb'^/12 + 2c/3. 

In what follows our fiuid has density 1, whereas V has density (J,0<(J<1((J = 1 and 
(7 = correspond to the uninteresting cases excluded above). In this setting Archimedes 
had to deal with three conditions: First, the floating condition V = aV (this condition 
is also known as Archimedes' principle, ibid. Sect, i), which, because of the above values 
of V and V, can be written 

a' = ay/a. (3) 
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Second, the equilibrium condition, which says that B — B' must be perpendicular to the 
plane £. Since we know (6, 0, —1) _L the equilibrium condition comes down to one of 



6 = or 56Vl2 + 2(c-a)/3 + 1/2 = 0. 



(4) 



Finally, the condition for the archimedean case, which reads 

a > —b^li + c. 



(5) 



In our formula language these conditions are not difficult to handle. However, one should 
be aware of the fact that Archimedes had no formulas at all but could only work with 
geometric propositions, which were enunciated in a rhetorical manner. In the case 6 = 0, 
([2]) yields c = a', and ([3]) gives c = a-^o". Since < a < 1, we see that c is positive and 
that ([5]) holds automatically. If 6 7^ 0, (E]) and (E]) yield 



and, thus, 



becomes 



c = a\fo — 6^/4, 



6^ = y(l-V^)-2. 



(6) 



(7) 



This can only hold ifa>3/4, soa<3/4 necessarily requires 6 = 0. Suppose now that 
6 < is a solution of (JTj). As ([5]) must be true, we obtain —b^/a < a(l — -y/o") + 6^/4 from 
(E]), and (I7j) transforms this inequality into 

-6v^ < 5678 + 3/4. 

This is true whenever a < 15/8. In the case a > 15/8, it is the same as saying that one 
of , 

-4v^ 



-^Ja /16a 6 

> 6 > — — + \ or 6 < 

-5 V255 -5 



16a 6 
"25" ~ 5 



holds. If we use ([7]) in the shape y/a = 1 — 3(6^ + 2)/(8a), we can read these inequalities 
as conditions for ^/a, namely. 



13 
cr > — 
" 25 



10a 5 V 25 



10a 



or 



^ 13 _ 3 _ 6 
^ - 25 ^ 10a ~ 5 



4 
25 



10a 



Summarizing we may say that a right hand equilibrium position in the archimedean 
case is possible only if a < 15/8 or if a satisfies one of the inequalities of ([H]). In these 
cases either 6 = or 6 < can be read from ([7]), whereas c is given by These formulas 
involve only rational expressions in a and a or square roots of such ones. Accordingly, 6 
and c can be constructed by means of ruler and compass if a and a are given. Archimedes 
fully described these constructions. The equilibria defined in this way can be classified 
as stable or unstable by means of the potential function, see Section 4. 

If we consider V (1 {z > bx + c} as the submerged part of V, our right hand position 
turns into a left hand one. Since V — V is the volume of the submerged part, the floating 
condition now reads a' = a\/a* for a* = 1 — a. In order to obtain an archimedean 
equilibrium, a must be replaced by a* in (^^| and ([H]). It is well known that the nature 
of the equilibrium (stable or unstable) remains the same, see fBj. 
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2. Equilibria in the non-archimedean case 



The case when £ intersects the basis circle {x^ + y'^ < a, z = a] oi V in more than one 
point was not considered by Archimedes. Henceforth it will be called the non-archimedean 
case. Note that neither case excludes the other, i.e., for a given paraboloid segment V 
with given density a archimedean equilibria may occur together with non-archimedean 
ones. One can also characterize the non-archimedean case by the condition 

a = hX + c with - i/a < X < ^/a. (9) 

The quantity X is closely connected with the size of the submerged part of the basis 
circle of P, so it has a natural meaning. Moreover, for most formulas it is advantageous 
to use X instead of c, so the reader should get accustomed to the fact that, by ©, the 
plane £ depends on h and X henceforth, i.e., £ = {z = a + b{x — X)}. Since formulas 
are slightly simpler for a left hand position, we assume that 

V' = Pn{z>bx + c} 

is the submerged part of V. The volume of V' equals Vi — V2, where Vi and V2 are the 
volumes of 

Vi=Vn{x>X,z<a} and V2 = V n {x > X, z < bx + c}, 

respectively. Here Vi and V2 are no longer paraboloid segments. We call Vi a right and 
V2 an oblique paraboloid sector. It turns out that the volumes of these sectors are closely 
related. 

In the case of the right sector Vi, a cross-section Pi fl {x = Xq} is a parabolic segment 
with area 4(a — x'^)'^^'^ /3. Hence 

4 f^, .,0/. a"^ fir X\ 2X^-5aX r- , , 
V, = -J^ (a-.r^<(. = -(^--arcta„^j+ ■/! (10) 

with X as in dH]) and 

A = a-X^ >0. (11) 

As a function of a and X, V\ will be denoted by V[a, X\ In order to obtain the volume 
V2, we consider another right paraboloid sector, namely, 

V[ = {x^ + y'^ < z, X > X\ z < a'}, 

with a' = 6^/4 + c defined as in ([2]) and 

X' = X~b/2. (12) 

By (Q, a' may be written as 

a' = 674 -bX + a. (13) 

The identity a' = X'"^ + A gives —\/a/ < X' < \/a/ (so V[ is really a right paraboloid 
sector). It is also useful to note 

A = a-X^ = a'- X'\ (14) 
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Clearly, V[ has the volume V{a',X') in the above sense. It is not hard to check that the 
affine mapping 

{x, y, z) ^ (x, y,z + hx) + (6/2, 0, li^/A) (15) 

induces a bijection between the right sector V[ and the oblique sector V2- Since the 
linear part of this mapping has determinant 1, V'l and V2 have the same volume 



a'2 / vr 



X' \ 2X'3 - 5a'X' 



V2 = V{a',X') = — [-- arctan--=j + VA, (16) 



where we have used f |T4l) . Now the floating condition (i.e., the analogue of ([3])) reads 

aV = Vi- (17) 

with V = a^7r/2 as in ([T]) and Vi, V2 as in (ITOl) . (fT6ll . respectively. 

Next we consider the center of gravity i?i = {xi,0,zi) of Pi. Its coordinates are 
defined by the moments 

XiV^i = - / x(a - x^f'dx, z^Vr = - / ^ (a - x^)^/^^^. (18) 
3 3 5 

This is easy to see, since the above cross-section Pi fl {x = xq} has the center of gravity 
(xo,0, (3a + 2xo)/5). Now one can verify (by differentiation of (ITSl) with respect to X, 
say) that these moments take the values 

xiV^i = ^A''\ z,V, = ^^1 + ^A^/l (19) 
15 6 45 

Obviously, these identities hold, mutatis mutandis, for the right paraboloid sector V[ as 
well. In particular, V[ has the center of gravity B[ = {x[, 0, z[), and the analogue of (fT^ 
reads 

x[V2 = ^A^l\ z[V2 = ^\/, + ^A'/' (20) 
15 6 45 

(recall A = a' - X'\ by i^). 

The affine mapping of (11 5p transforms into the center of gravity B2 = (x2, 0, Z2) 
of the oblique paraboloid sector V2, so we have 

X2 = x'l + 6/2, ^2 = -z'l + bx[ + 6^/4. 

Together with (1201) . ffT^ and ffT^ this gives the following formulas for the moments X2V2, 
Z2V2: 



T. 4 /562 2bX 2a\ f AX 2b\ 



(21) 



Let (x', 0, z') be the center of gravity of V. Obviously, our moments satisfy 

x\Vi - V2) = xiVi - X2V2, z\Vi - V2) = ziVi - Z2V2. (22) 

We are now in a position to enunciate the analogue of (jl]) in the non-archimedean case. 
As B — B' must be perpendicular to we obtain 

2a/3 — z' = x'/6. 
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Since the volume Vi — V2 of V' is positive, this is equivalent to 



(2a/3 - z') {Vi - V2) - (x'/b) {Vi - V2) = 0. (23) 

Here we insert the right hand sides of (12^ for x'{Vi — V2) and z'{Vi — V2). Moreover, we 
use (fTOjl and (|2Ti) . Then a short calculation shows that the left side of (1231) comes down 
to 

^ /562 26X 1\ 2b , 



^12 3 27 9 
(here stands for equilibrium ). We recall (ITTI) and summarize our results in 

Theorem 1 In the non-archimedean case a left hand equilibrium position is character- 
ized by the conditions F = and E = 0, where 

F = Vi-V2-aV 

and E is given by (p^ . The quantities V , Vi, V2 and A are as in ([1]), (fTOj) . (fT6|) and 



Remarkably, the equilibrium condition E = involves only the volume V2, whereas the 
floating condition F = also involves Vi. Further, the expression 

/ = 55712 - 2bX/3 + 1/2 (25) 

occurring in E is, by virtue of ([H]), identical with the expression 56^/12 + 2(c — a)/3 + 1/2 
from the equilibrium condition (jl]) in the archimedean case. It seems, however, that this 
formal analogy has no influence on the rather different properties of both cases. 

For a right hand position of V and density a, the floating condition reads Vi — V2 — 
a*V = with a* = l — (see the Archimedean case). The equilibrium condition remains 
unchanged. 



3. Finding "all" solutions in the non-archimedean 
case 

As above, suppose a and a are given. Then we know that V can take at most two 
archimedean right hand equilibrium positions, corresponding to 6 = or 6 = — (8a(l — 
^/o")/3 — 2)^/^, see (HI), ([7]). In the non-archimedean situation the determination of the 
number of possible equilibria in a mathematically rigorous way seems to be much more 
difficult, and even a general upper bound for this number is out of reach for us. In fact, 
we do not know how to bound the number of zeros of the (non-algebraic) map 

] - 00, 0[ X ] - v^, v^[^ M2 . ^) ^ ^26) 

If, however, we assume that X (instead of a) is given together with a, we can deter- 
mine the exact number of values b satisfying E = quite well (see Theorem [2]). Each 
of these values b gives, because of F = 0, exactly one density a. In this way we obtain 
reliable diagrams connecting X and a. For this reason a reliable (though not rigorous) 
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answer to the question about the number of zeros of (|26ll seems to be possible. As we 
remarked aheady, X can stand for the size of the submerged part of the basis circle of 
V] hence replacing a by X is not quite unnatural but has a certain value of its own. 
Recall that E = fV2 + 2bA^/^/9 with / as in Put 

~ E 1 /tt X'\ 2X'^ - 5a' X' r- 2b A^/^ 
E = — — - = arctan —== H v A H — — , 

with a', X', A as in (fT3|) . f|T2l) . (fTTj) . respectively. The advantage of E lies in the fact 
that its derivative with respect to 6 is a rational function of b. Indeed, 

dE ^ ^5/2 



db 108a'3/2' 

where P is the polynomial 

P = 6Xb^ + (-10a + 21) 6^ - 36X 6 + 12a + 18, (27) 

which is cubic in b. We further note that / has zeros bi < 62 < only if X < -^15/8. 
In this case these zeros read 

4X 1 , — ^ , AX 



61 = — - -V16X2-30, 62 = — + -V16X2 - 30. (28) 
5 5 5 5 

With these tools at hand, we are able to enunciate the main result of this section, which 
describes the (negative) solutions b of E = 0. 

Theorem 2 Let a > and X be given, —\fa < X < ^a. 

(a) If X < — a/15/8, P has exactly two negative zeros 61 < 62- Moreover, bi < bi, 
and E = has exactly two solutions, one in ]bi,bi[, the other in ]b2,0[, where bi, 62 o'^e 
as in 



(b) In the case —^/iKJS < X < the equation E = has solutions only if P has two 
negative zeros bi < 62- If E{b2) = 0, then 62 is the only solution of E = 0. If E{b2) < 0, 
then E = has exactly two solutions, which lie in the intervals ]6i,&2[ (^nd ]&2,0[. // 
£'(62) > 0, E = has no solution. 

(c) In the case X = the equation E = has no solution for a < 21/10. If a > 21/10, 
E = has exactly one solution, which lies in ]&i,0[, where bi is the negative zero of P. 

(d) In the case X > 0, P has exactly one negative zero bi and E = has exactly one 
solution, which lies in ]&i,0[. 

Our Proof of Theorem [2] requires the knowledge of the limits of E for b tending to or 
to —00. We note 

Proposition 1 In the above setting, 

lim£; = — >0 and lira E = ^^A^^^ . 

b^O 2 45 
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Proof. The value of the first hmit is clear from (|2^ and the fact that V2 tends to Vi for 
6 — *• 0. As to the second one, recall that V2 has the center of gravity (x2, 0,2:2), where 
X2, Z2 satisfy (j2Ti) . Now the second identity of (|2Ti) can be written 

Since V2 ^ for b —00 and < Z2 < a, we conclude that Z2V2 tends to zero and E to 
the value in question. ■ 

Proof of Theoreml^ Note a' > for all X in question, so only the zeros bi < &2 of / (as 
given in fl28l) ) can be poles of E. 

This situation occurs in case (a), which we discuss first. The definition of E shows 

lim E = 0, lim E = —00. 

b— >— 00 6— >bi,fe<fei 

By Proposition [H E must be positive for all 6 -C 0; hence it takes a positive maximum in 
] — 00, bi[. This requires that dE/db and, thus, P, has a zero 61 < fei, for which the said 
maximum is taken. By Descartes' rule (see [7], p. 310), P has exactly one positive zero, 
so it must have another zero 62 < 0. We may assume 61 < 62- In the case X = — a/15/8 
one verifies 62 = ^1 = ^2- If < — a/15/8, we have 61 < 62 and 

lim E = 00 = lim i?; 

6— >bi,fe>bi b^b2,b<b2 

from this we conclude that ii^ takes a minimum in ]6i,62[, more precisely, for 6 = 62 G 
62 [• In this way we know the intervals where E is strictly monotonous. We see, first, 
that E has exactly one zero < 61, which lies in ]bi, bi[. Since 

lim E = — CX3, limii^ > 

b^b2,b>b2 b^O 

(recall Proposition [1]), we see, second, that E has exactly one zero in ]&2,0[. As E is 
positive in ]6i,62[, our list of (negative) zeros of E is complete. 

In the remaining cases poles of E occur no longer. In case (b) the polynomial P also 
has exactly one positive zero. If P has no zeros < 0, dE/db is positive in ] — 00, 0[, so E 
is strictly increasing. Since both limits of Proposition [T] are positive, E does not vanish 
in ] — 00, 0[. This assertion remains true if P has a double zero < 0. Suppose, therefore, 
P{bi) = = -P(&2) for 61 < 62 < 0. It is not hard to see that E must have a positive 
maximum at bi and a minimum at 62- Monotonicity arguments and consideration of the 
sign of -^(62) yield the number and the location of the solutions of E = just as indicated 
in the theorem. 

In case (c), P becomes a quadratic polynomial, which is positive for all 6 < as long 
as a > 21/10. In this case E is strictly increasing in ] — 00, 0[, and since it tends to zero 
for b —00, it must be positive throughout. If a < 21/10, P has exactly one zero bi in 
] — 00, 0[. Now E is decreasing for b < bi but increasing for 6 > 61, so it takes a minimum 
for b = bi. The limit of E for b —00 being 0, we conclude E{bi) < 0. These arguments 
imply that E has exactly one zero < 0, which lies in 0[. 
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In case (d), Descartes' rule, when applied to the polynomial P{—b), shows that P has 
exactly one zero bi < 0. Since E < for all 6 -C and limb_»_oo E = 0, the function E 
has a negative minimum at hi. Therefore, E vanishes for some b g]6i,0[, and this b is 
the only zero of E for reasons of monotonicity. ■ 

Example 1. We apply Theorem [2] in the case considered by Rorres [8]: He took a = 
0.51 and a "base angle" ip = 74.33. This value of ip corresponds to a = (tan(y9)^/4 ^ 
3.17690918 in our setting. From Theorem[2]we know that, for each X e] — y/a, — a/15/8 ], 
the equation E = has exactly two solutions b e] — oo,0[. Hence we expect that 
E = describes a curve (X, b) with two branches for all X in this range (note —^/a ^ 
-1.782, -^15/8 ^ -1.369). Next our theorem suggests testing whether E{b2) < for a 
sufficiently large number of (equidistant) values of X in ] — a/15/8, 0[. Since this is true 
in all cases, the curve (X, b) will also have two branches for X in this interval. Finally, 
the theorem implies that (X, b) consists of one branch for X G [0, ^/a[. Using the floating 
condition F = of Theorem [1] in the shape a = (Vi — V2)/V, each pair (X, b) produces 
a pair (X, a). One expects that each branch of (X, 6) produces a branch of the curve 
(X, a) in this way. 

In order to obtain Diagram 2, one has to investigate all values X between —1.78 and 
1.78 in steps of 1/100. Numerical values for the corresponding solutions b of E = can be 
found by standard methods as Newton's algorithm. The diagram displays the resulting 
points (X, a). Our steps are small enough to produce the picture of a nearly continuous 
(and smooth) curve — in contrast to the situation of Example 2 below. The horizontal 
line a = 0.51 intersects the curve in one point of the upper branch and in at most two 
points of the lower one. Hence we expect at most three left hand equilibrium positions 
to be detected in this way. 




Inspecting our computation more closely we find the pairs 

(X,a) ^ (-1.04,0.50997999), (-1.03,0.51000418), (-1.02,0.50999785) 



10 



on the lower branch. Accordingly, this branch should contain two equilibria with a = 
0.51. With these starting values and two more on the upper branch it is no more difficult 
to find 

{X,b) ^ (-1.03304236, -1.12424322), (-1.02105684, -1.13986072) 
and (-0.12106085,-12.68795681) 

(these numerical values of a, X and b satisfy E = and cr = 0.51 up to an error < 10^^). 
The corresponding approximate "tilt angles" (in Rorres's sense) are 131.653°, 131.260° 
and 94.506°, respectively. For right hand equilibria we have to replace a by a* = 1 — a, 
which means that we intersect our curve with the horizontal line a = 0.49. In this way 
we detect two more solutions 

(X, b) ^ (-1.46372405, -0.69920557), (-0.74316119, -1.52773443), 

which correspond to approximate tilt angles of 34.961° and 56.793°. Altogether, we have 
found five non-archimedean equilibria in this case; Rorres has only four, since he considers 
those corresponding to angles of 131.653° and 131.260° as only one solution with an angle 
of 131.5°. Rorres' situation can be established by a slight modification of a. Indeed, 
a ^ 0.51000554 melts those two solutions into (X, 6) ^ (-1.02702703,-1.13205421), 
the corresponding tilt angle being ^ 131.456°. However, Rorres writes a = 0.510 and 
ip = 74.330, which suggests interpreting these figures as exact values. 

Remarks. 1. The reader need not worry about the fact that a change of the last digit 
of our values for X or 6 may give slightly better results in combination with the above 
numerical approximation 3.17690918 of a. Our results have been obtained by rounding 
off the higher digits of substantially better values — not only of X and b but also of a; 
accordingly, X and b may appear not sufficiently precise when they are combined with 
this approximation of a. 

2. Computations suggest that the number of five non-archimedean equilibria (as in 
the example) could represent the maximum. 

In this example each value X g] — ^/a, ^/a [ defines at least one point (X, b) that satisfies 
the equilibrium condition E = 0. This is no longer true if a < 3. In this case an 
investigation of the discriminant of the polynomial P of (127|) exhibits values a > and 
X < for which E is monotonically increasing; by Proposition [T], this means that E = 
has no solution 6 g] — oo, 0[. We summarize these observations in 

Proposition 2 Let a > and put ai = (-213 + 198v/Tl)/250 1.7748;. Then the 
equilibrium condition E = has no solution 6 g] — oo, 0[ if 

(a) a < ai and X g] — y/a, 0[, 

(b) a G]ai, 21/10] and X G [Xi,0[ or 

(c) aG]21/10,3] andX e [Xi,X2], 

where Xi = -(7 + ^5/27)^/^ X2 = -(7-75/27)^/2 ^^th = -llaV54 + 5a/9 + 13/24 
and 5 = {3-a){a + 6f. 

Proof. We briefiy sketch the main arguments. Suppose X < 0. Then P has a zero > 0. 
Further, E is monotonically increasing (on ] — 00, 0[) if, and only if, P has either no zero 
< or a double zero < 0. This is the same as saying that the discriminant D of P is 
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< 0. But D can be written as a quadratic polynomial Di m Y = X^. The complex 
zeros of D^ are Yi = -f + 75/27 and 1^2 = 7 - 75/27. These zeros are real only if a < 3, 
which means that only in this case D < is possible. Assuming a < 3 henceforth, one 
sees that Yi is always positive, whereas Y2 is positive just if a > 21/10. This gives the 
intervals for X in the cases (b) and (c). However, one has also to observe that Yi < a if, 
and only if, a> ai (whereas Yi < 15/8 is always true). ■ 

Example 2. Consider a = 5/2, which falls under case (c) of the proposition. Accordingly, 
the curve (X, 6) is empty for X e [Xi.Xi], where Xi ^ -1.143 and X2 ^ -0.0917. 
Since a/5/2 1.581 and a/15/8 ~ 1.369, we know from Theorem [2] that it should have 
two branches for —1.58 < X < —1.37. By means of this theorem one verifies that these 
two branches should extend as far as X = —1.29. Furthermore, two branches are to be 
expected for —0.08 < X < 0, but only one branch for X > 0. 

The diagram displays points of the curve (X, a) with —1.58 < X < 1.58, again in steps 
of 1/100. A comparison with Example 1 reveals some marked differences. First, the two 
branches for —1.58 < X < —1.29 are connected, and the same holds for —0.08 < X < 0. 
Second, the diagram suggests that there is exactly one non-archimedean solution for 
0.415 < a < 0.585 (observe that possible right hand equilibria with —1/2 < a < 0.585 
correspond to left hand ones with 0.415 < a < 1/2; for a = 1/2 see Section 6). The 
density of points is rather low when the tangent of the curve is nearly vertical, so a 
"continuous" picture requires much smaller steps for values of X in this region (say steps 
of 1/10000 instead of 1/100). 

tI 



I 

-1.58 



a = 0.585 



a 



0.415 



Diagram 3 




The considerable differences between our examples suggest that a global theory of the 
non-archimedean case (such as a theorem about the number of equilibria for a given pair 
(a, a)) may be a difficult matter. 



12 



4. Classification of equilibria in the archimedean case 



The classification of tlie above equilibria requires considering the potential energy of a 
certain position of the paraboloid segment V. As in Section 1, we start with a right hand 
position in the archimedean case. Recall that the center of gravity of the submerged part 
V is B' = {x', 0, z') with x' = 6/2 and z' = + 2c/3. We need the moments x'V 

and z'V , where V = a''^n/2 is the volume of V (with a' as in ([2])). In what follows we 
use the abbreviation 

p = VipTi. (29) 

Further, we work with a Hesse normal form of the plane S in order to describe the 
distance of a point (x, y, z) from £, namely 

^ = {(^_6a;-c)//? = 0}. 

Therefore, —{z' — hx' — c)V' / (3 can be considered as the potential of the buoyancy of V. 
In the same way the potential of the weight of V is given by (2a/3 — c)aV/(3 since V 
has the volume V and the center of gravity B = (0,0,2a/3), see ([1]). Accordingly, our 
potential function has the shape 

U = U{c, b) = ((2a/3 - c)aV + aV'/3)/(3 

because —z' + bx' + c = a'/3. In this section we write 

Fo = V'-aV and Eo = fV', 

where / = 56^/12 + 2(c — a)/3 + 1/2, see @. The floating and equilibrium conditions 
read, thus, Fq = and bE^ = 0, respectively. One easily verifies 

dU Fq dU b ( (2a \ ^ \ 

II _ c Fo + . (30) 



do (3' db f3^ Vv 3 

Equilibria should be the same as stationary points of the potential function. An inspection 
of the derivatives of U as given in ( l30l) proves this. The second derivatives of U can be 
written 

d^U 2V' d'^U b fV Fn 



dc^ a' 13' dcdb (3 \a' /52 



(31) 



d^U b^V f5b^ c+l\ l-2bWf2a ^. 

+ + — - c Fo + Fo . (32) 



and 

db^ a'f3^ V 8 ' 2 ; ' /35 3 

In the equilibrium case with 6 = we also have Fq = 0. Therefore, (l3Tll and ( |32l) show 
that the Hesse matrix of U has the shape 

2V'/a' 
Fo 

Recalling a' = > we see that this matrix is positive definite if, and only if, Fq > 0, 
which is the same as saying a < 3/(4(1 — -y/a)). So in this case the stationary point of U is 
a minimum and our right hand equilibrium is stable; it becomes unstable (more precisely, 
a saddle point) for a > 3/(4(1 — -^/o")), whereas the case a = 3/(4(1 — -^/o")) cannot be 
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classified in this way. Using higher derivatives one can show that this equihbrium is also 
stable (as Archimedes did without this device, see [1], Lib. II, Sect. iv). Similarly, the 
case b < 0, Fq = Eq = 0, gives the Hesse matrix 

_r / 2 b \ 

a'P V b (5678 + (c+l)6V2)//32 J " 

Again, a' = ay/a > 0. Since the determinant of this matrix equals V'^b"^ / {a' P^) > 0, we 
see that equilibria with 6 7^ are stable. 

Remarks. 1. In our model no maximum of the potential function is to be expected, 
since the potential energy will always grow if one moves V "upwards" , i.e., more to the 
right hand side and, simultaneously, in a direction perpendicular to S. This observation 
applies to left hand positions in an analogous way, for instance, to the non-archimedean 
positions of the next section. 

2. The left hand position that corresponds to our right hand one has the same 
potential, see [U] and the end of Section 1. This can be checked directly if one observes 
that a* = 1 — a plays the role of a and 

B" = B + ^^iB-B') 

that of B; observe, further, that the weight and buoyancy potentials change their signs. 



5. Classification of equilibria in the non-archimedean 

case 

As in Section 2, we consider a left hand position of V in the non-archimedean case and 
adopt the corresponding notations. In order to define the potential function U, we use 
the same Hesse normal form of S as as in the foregoing section. Then the potential of the 
weight of V remains unchanged up to the sign, i.e., it equals —{2a/3 — c)aV/P, with V as 
in ([T]) and (3 as in 0291) . The sign change is due to the transition from a right hand position 
to a left hand one. Similarly, the potential of the buoyancy is {z' — bx' — c)(yi — V2)/ [3. 
Hence we have 

U=-^{{z'-bx'-c) {y,-V,)-{^-^-<^ 

In the spirit of Section 2, we consider ?7 as a function of X and b instead of c and 6, so 
we write c = a — bX. By means of formulas (l22l) and (I2T1) we obtain 



U = UiX, = ^ ((I - bx) {aV - V,) + ^V, - ^) (33) 
(recall (fT3l) ). Now the analogue of ( j30l) reads 

au 6F au 1 /,,eUx^'-^)f) (34) 



dX (3' db (3^ \ V 3 
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with E and F as in Theorem [H We also note the analogues of (|3Tll and ( |32l) , namely, 

where X' = X - 6/2 (see ([I2D) and 
d^U 1 



-2X6^ + (4a - 7)6^ + 14X6^ - 66 + 12X)\/2 + + ^i, (36) 



with 



Fi = 4;-(-2a62 -9X6 + a)F and 

El = -^(46^-4X6^+(13-8a)6^-28X62 + (4a + 6)6-12X)E. (37) 

Of course, the correctness of these formulas is easy to check by the aid of a computer 
algebra system, say. However, it takes some effort to find them and we think, therefore, 
that it is justified to render them here. In the case of an equilibrium we have Fi = Ei = 0, 
so ( 1371) disappears and ( 136|) looks fairly simple then. This obviously happens for the 
second item of 0351) . too. Further, formulas 0341) show that there is hardly a simpler 
characterization of stationary points of U than our floating and equilibrium conditions. 

Example 3. We return to Example 1. Our formulas (l35l) and (|36l) quickly give the Hesse 
matrix of the equilibrium positions we described there. From altogether three left hand 
equilibria (with a = 0.51) the first and the third one (with tilt angles of about 131.653° 
and 94.506°) are stable, since the Hesse matrices have the pairs of approximate eigenval- 
ues (0.00101514, 6.83907084) and (0.00001567, 7.50021176) in these cases. The second 
left hand equilibrium (131.260°) is unstable (a saddle point), the respective eigenvalues 
being approximately (—0.00098808,6.78938522). In the same way the first of the right 
hand equilibria (34.961°) is stable and the second one (56.793°) a saddle point. The case 
when the first and the second equilibrium melt into one (with an angle of about 131.456°) 
can be settled by the aid of higher derivatives. For this purpose we write (Xq, 6o) for the 
corresponding value f» (-1.02702703,-1.13205421) of (X,6). We use the substitution 
X = F + A6, where 

d^U/dXdb 

On differentiating U with respect to Y and 6, one sees that the second derivatives vanish 
for the respective point with the exception of d'^U/dY'^, whereas d^U/db^ takes a value 
c ~ 0.20378903. Hence the behaviour of [/ in a neighbourhood of (Xq, 6o) is like that 
of cZ^ + 0(y2 -\- ^YlZ"^ + Z^) for Z, Y close to zero. This expression, however, becomes 
negative for Y = and Z < 0, \Z\ small. Therefore, (Xq, 6o) defines an unstable position. 

Of course, we can follow the same line when we investigate the stability properties 
of those points on the curve (X, a) that are rendered in Diagram 2. We find that all 
points on the upper branch describe stable equilibria. As to the lower branch, those with 
X < —1.03 belong to stable equilibria and the remaining ones (with X > —1.02) to 
saddle points. It seems, thus, that (Xq, ctq) with (Tq ~ 0.51000554 (which comes from the 
above pair (Xo,6o)) forms a sort of limit point for the stability of left hand equilibria. 
Additional computations (with smaller steps) confirm this observation. 
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Again, the behaviour of Example 2 is different. All points of Diagram 3 with X < 
— 1 (i.e., those of the left component of the curve) give stable equilibria, whereas the 
remaining ones belong to saddle points. 



6. The horizontal case 

The case when the axis of the paraboloid segment is horizontal was not treated so far. 
It will henceforth be called the horizontal case. Our attempts to include this case in the 
above discussion gave rise to problems with differentiability, so it seems justified not to 
do so. In our setting the horizontal case can be characterized hy £ = {x = X}, where X 
satisfies —\/a < X < ^a, see ([9]). One easily checks that the equilibrium condition takes 
the simple shape (2a/3 - zi)Vi = 0. By this is equivalent to 4X^^/^/45 = and, 
therefore, to X = 0. In this case Vi = V/2 holds for the respective volumes, so the floating 
condition Vi = aV requires a = 1/2. Accordingly, only this rather obvious equilibrium 
position is possible in the horizontal case. The classification of this equilibrium, however, 
is less obvious. 

To this end we use the potential Uq of this position; it takes the value 4a^/^/15, as is 
readily seen. Moreover, we consider neighbouring left hand positions, i. e., pairs {X,b) 
with |X| small and —b large. It is advantageous to work with X' = X — 6/2 ^ instead 
of b (see ( !T^ ). Therefore, the potential function U of flH^ reads U = U{X,X') now. 
Then we insert the series 



X X X^ X^ X' TT VA ^3/2 
arctan —= = —= -—7- H -— ■ — . . . , arctan —= = h 



and 



y/A y/A 3/13/2 5^5/2 ^ 2 X' 3X'3 



1 X X2/2-1/16 



/5 v/4(X' -Xy + 1 2X' 2X'2 X'3 

together with 



1 1 5X2 35jj^4 

+ + 



a5/2 2a7/2 8a9/2 
into flHHj) and obtain, in a straightforward (though laborious) way, 

2^3/2 45/2 /AflVS fl5/2\ 

f/(X, X') = Uo + ^X2 + ^XY + ^]y^ + 0((|X| + Yf), 

^ ' ^ " 3 15 V 105 30 y vvi I ; 

where Y = 1/X' is positive and close to 0. The quadratic form in X, Y on the right 
hand side is positive definite for a > 35/12 and indefinite for a < 35/12, thus indicating 
stability and instability, respectively. To settle the case a = 35/12, we need more terms 
of this expansion and the substitution Y = Y' — 12X/7. This gives 



u = Uo- !^^x^ + o{x'\Y'\ + + 1^|5^ 

gggg ^ V I I ^ ^\ \ ) 



(where Uq takes the value 245^105/648). So we can choose, for each small value of |X| 
a number Y' with \Y'\ < \X\ such that U < Uq. Hence this equilibrium is unstable. 



16 



References 

[1] Archimedes, De corporibus fluitantibus, in: Opera omnia, vol. 2, 2nd ed. by J. L. 
Heiberg, Teubner, Leipzig 1913; 317-413 

[2] Archimedes, On floating bodies, translation by T. L. Heath, in: Great Books of the 
Western World, vol. 11, Encyclopedia Britannica, Inc., Chicago 1952; 538-560 

[3] Archimedes, Uber schwimmende Korper, Ubersetzung von A. Czwalina, in: Werke, 
Wissenschaftliche Buchgesellschaft, Darmstadt, 1972; 283-346 

[4] A. Czwalina, Die Mcchanik des schwimmenden Korpers, Akademische Verlagsge- 

scllschaft, Leipzig 1956 

[5] E. J. Dijksterhuis, Archimedes, Munksgaard, Copenhagen 1956 

[6] E. N. Gilbert, How things float, Amer. Math. Monthly 98 (1991), 201-216 

[7] A. 1. Kostrikin, Introduction to Algebra, Springer, New York 1982 

[8] C. Rorres, Completing book II of Archimedes's On Floating Bodies, Math. Intell. 
26 (2004), 32-42 

[9] S. Stein, Archimedes and his floating paraboloids, in: Mathematical Adventures, ed. 
by D. F. Hayes and T. Shubin, Math. Assoc. America, Washington 2004; 219-231 



Kurt Girstmair and Gerhard Kirchner 
Institut fiir Mathematik 
Universitat Innsbruck 
Technikerstr. 13/7 
A-6020 Innsbruck, Austria 

Kur t . Girstmair @uibk . ac. at , Gehard . Kirchner @uibk . ac. at 



17 



